Analysis of the elementary excitations in high- Tc cuprates: 
explanation of the new energy scale observed by ARPES 
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Using the Hubbard Hamiltonian we analyze the energy- and momentum- dependence of the ele- 
mentary excitations in high-Tc superconductors resulting from the coupling to spin fluctuations. 
As a result of the energy dependence of the self-energy ^(k,oo), characteristic features occur in the 
spectral density explaining the 'kink' in recent ARPES experiments. We present results for the 
spectral density A(k,a;) resulting from the crossover from Im S(k,a;) oc a; to Im S(k,a;) oc cj^, for 
the feedback of superconductivity on the excitations, and for the superconducting order parameter 
A(k, (jj). These results relate also to inelastic neutron scattering and tunneling experiments and shed 
important light on the essential ingredients a theory of the elementary excitations in the cuprates 
must contain. 
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For understanding the high-Tc cuprates their ele- 
mentary excitations are of central significance. Angle- 
resolved photoemission spectroscopy (ARPES) is a pow- 
erful tool for studying the elementary excitations in 
high-Tc superconductors because the spectral density 
contains all information about self-energy effects. Due 
to an improved angular resolution (momentum dis- 
tribution curve (MDC) and energy distribution curve 
(EDC)), data became available which provide new in- 
sight on the momentum and frequency dependence of 
the self-energy i;(k, cj). In particular, a 'kink' feature 
at fio; ~ 50 ± 15 meV has been observed in hole-doped 
cuprates like Bi2Sr2CaCu208, Pb-doped Bi2Sr2Cu06, 
and La2-a;Sra;Cu04 Experiments by Shen et at., 

[||,|| observe the kink feature in all directions in the first 
Brillouin Zone (BZ). It exists in both the normal and su- 
perconducting states. On the other hand, Kaminski et 
al., [|j discuss the break only along the (0,0) (tt, tt) 
direction occuring when one goes from the normal to the 
superconducting state. Therefore, they did not analyze 
the feature observed by the other group [^,||. However, it 
is quite interesting that a close analysis of data of Kamin- 
ski et al. 1^ in the normal state reveals the same changes 
of the Fermi velocity, vp^ as noted by Shen et al. [^,^. 
Thus, there seems to exist a 'new' energy scale in hole- 
doped cuprates. Remarkably, the electron-doped coun- 
terparts {e.g. Ne2-xCeajCu04) do not show a 'kink' [||. 
So far, interpretations are given in terms of the presence 
of a strong electron-phonon interaction [||,||, stripe for- 
mation 0, or coupling to a resonating mode [HJ|]. It is 
interesting that the experiments also observe a change in 
the dispersion of the elementary excitations going from 
the normal to the superconducting state We will 

show that this results from the feedback effect of super- 
conductivity on the elementary excitations. 

In this letter we present a study of the spectral den- 
sity A(k, a;) of the elementary excitations using an elec- 
tronic theory based on Cooper-pairing due to an ex- 
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FIG. 1. Calculated self-energy effects in the spectral den- 
sity of the quasiparticles in hole-doped superconductors in the 
normal state at T = lOOK and along the (0, 0) (tt, 0) direc- 
tion. The dashed line at = denotes the unrenormalized 
chemical potential. In the inset the spectral density along the 
(0,0) (tTjTt) direction is shown. In both cases at energies 
about /icj 65 =b 15meV a 'kink' occurs, since the velocity of 
the quasiparticles changes. The results are in good agreement 
with experiments (see, for example Fig. 3 in Ref. |Q or Fig. 4b 
in Ref. [^). Note, the width of the spectral density peak is 
for (0, 0) (tt, 0) twice the one for (0, 0) (tt, tt). 

change of antiferromagnetic spin fluctuations. In par- 
ticular, we show that the 'kink' in the spectral density 
can be naturally explained from the interaction of the 
quasiparticles (holes) with spin fluctuations. In agree- 
ment with recent experiments we will demonstrate that 
the 'kink' feature is present in both the normal and su- 
perconducting state [^,^-^. Thus, we are able to ex- 
plain recent ARPES experiments which study in detail 
the spectral density and in particular the energy disper- 
sion C(;(k) = e(k) + I](k, cj). It is significant that the 
self-energy I](k, cj) resulting from the scattering of the 
quasiparticles on spin fluctuations can explain the main 
features observed. We argue that our results for the ele- 
mentary excitations suggest a crossover from Fermi liquid 
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FIG. 2. Positions of the peaks in the spectral den- 
sity A(k,co) versus k-k^ (energy dispersion) along the 
(0, 0) (tt, 0) direction of the BZ calculated within the 
FLEX approximation. This has to be compared with the po- 
sition of the peaks derived from the momentum distribution 
curve (MDC) for a hole-doped superconductor as measured 
in experiment. The curves show a 'kink' at energies about 
hou ^ 65 ± 15meV. The dashed line is a guide to the eyes. 
We find small changes due to superconductivity which almost 
coincide with the 'kink' position. Inset: change in the peak 
position in A(k, oo) in the superconducting state (T = 0.5Tc). 
The results are in fair agreement with ARPES data 

to a no n- Fermi liquid behavior. Furthermore, we obtain a 
picture consistent with inelastic neutron scattering (INS) 
and tunneling measurements 

The theoretical analysis is based on the spectral den- 
sity for the elementary excitations within Nambu 
space [|ll| which are given by, after continuation to the 
real cj-axis. 



A(k,cj) = 



S''(k,cj) 



(1) 



TT [c^-6k-S'(k,Cj)]" + p''(k,C^)]" 

Here, ek is a tight-binding energy dispersion on a square 
lattice and i;'(k, c^;) and i;''(k, cj) are the real and imag- 
inary part of the self-energy, respectively. We perform 
our calculations for the elementary excitations 



c^(k,T) = e(k) + S(k,cj(k,T),T) 



(2) 



The superconducting gap function ^(k, cj) is calculated 
self-consistently using the 2D one-band Hubbard Hamil- 
tonian for a Cu02-plane, which reads on a square lattice 



H 



(3) 



Here, cf^ creates an electron with spin a on site z, U 
denotes the on-site Coulomb interaction, and tij is the 
hopping integral. The imaginary part of the self-energy 
is given by Q 

f/2 



Im I](k, uj) 



coth 



2T J 



tanh 



2T 



^Imx(k-k^c^O S{\u;-co^\-e^^) 



(4) 



where Im xCq?^) is the imaginary part of the spin 
susceptibility within the random phase approximation. 
We determine the coupling of the quasi-particles to the 
spin fluctuations using an effective perturbation theory 
(FLEX) [l2-14| which we calculate directly on the real 
cj-axis. 

These equations are standard, however it is important 
to realize that due to the combined effects of Fermi sur- 
face topology and x(q = Q,^) at the antiferromagnetic 
wave- vector Qaf = (7r,7r), the k and uj dependence of 
I](k, a;) become very pronounced and change the disper- 
sion a;(k). It is known that the strong scattering of quasi- 
particles on antiferromagnetic spin fluctuations results 
in a non-Fermi liquid behavior of the quasiparticle self- 
energy for low-lying energy excitations, in particular, in 
Im ^ CO Clearly, it follows already from the 

Eq. (H) that the expected doping and momentum depen- 
dence resulting from the crossover from S oc o;^ to S oc cj, 
i.e. to a non- Fermi liquid behavior, can be reflected in 
cc;(k) and yl(k, cj). Simply speaking, the change in the 
cc;- dependence of the self-energy i;(k, cc;) changes the ve- 
locity of the elementary excitations. Thus, for a given 
k-vector, the MDC curve shows a 'kink' at some char- 
acteristic frequency controlled by Usf (cc;s/=spin fluctua- 
tion energy). Regarding the superconducting state the k- 
and cc;-dependence of the order parameter A(k, cj) is im- 
portant and yields the feedback of the superconducting 
state on the elementary excitations. 

This quite new structure in a;(k, T) which is present 
in both the normal and superconducting state is shown 
in the figures exemplarily for optimal doping and results 
from our calculations obtained by solving the above equa- 
tions self-consistently within a conserving approximation 
[|l3[. The full momentum and frequency dependence of 
the quantities is kept and no further parameter is intro- 
duced. 

In Fig. we present results for the frequency and mo- 
mentum dependence of the spectral density in the normal 
state exemplarily along the (0, 0) (tt, 0) direction cal- 
culated using the canonical parameters U = At^ and t = 
250 meV [0. The changes in the k-dependence of the 
peak in A(k, cj) reflect the characteristic features in the 
self-energy S(k, cj) or in the velocity of the quasipar- 
ticles. The kink occurs at energies about fio; ~ 65 ± 15 
meV for optimal doping (x = 0.15) and Tc ~ 65K. We 
also find that the 'kink' feature is present in all direc- 
tions in the BZ {uo ^ uosf + Vi?(^)k, k=k(A:,0)) and, in 
particular, along the diagonal (0, 0) (tt, tt) direction 
as shown in the inset of Fig.l. We get that the kink is 
similar pronounced in both directions. Moreover, we see 
from our calculations that this feature has only a weak 
temperature dependence over a wide temperature range. 
It changes only at very small temperatures which we will 
describe later. 

In Fig. 2 we show the positions of the peaks along 
(0, 0) (tt, 0) shown in Fig. 1 as a function of (k — ki?) 
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FIG. 3. Calculated frequency dependence of the quasipar- 
ticle self-energy at the node S(kn,a;), kn = (0.4, 0.4)7r. The 
solid curves correspond to the normal state at T = 2Tc, 
whereas the dashed curves refer to the superconducting state 
at T = 0.5Tc. At k = kn, where the superconducting gap 
vanishes, one clearly sees approximately at Tiix) — 65meV 
a crossover from Fermi liquid behavior (S oc o;^) to a 
non-Fermi liquid behavior (S oc uo) for low-energy frequen- 
cies as a function of temperature. We show in the inset 
the behavior of E(kn,a;) calculated at very low temperatures 
T = 0.003t 0.9K (dashed line). 

for different temperatures. We obtain only small changes 
due to superconductivity which almost coincide with the 
'kink' position. Remarkably, the deviation at k — k^? ^ 
0.05 is due to the frequency dependence of the self- 
energy and reflects the transition from Fermi-liquid to a 
non-Fermi liquid behavior along the route (0, 0) (tt, 0) 
in the Brillouin Zone (BZ) for both normal and super- 
conducting state. In the inset we show results for the 
difference in the peak positions for the normal and su- 
perconducting state along (0, 0) (tt, 0) direction in or- 
der to see also the feedback of superconductivity. Note, 
this disappears for k — kj? ^ 0.05 corresponding to ap- 
proximately 65meV. We get changes of about 10 meV 
while Kaminski et al, Q observe along (0,0) (tt, tt) 
direction a larger difference of about 20 meV. By inspect- 
ing Eq. (||) this is expected, since the feedback effect of 
superconductivity on x is larger for (7r,7r). 

In order to investigate the effect of the self-energy 
S(k, cj) on the dispersion a;(k, T) we show in Fig. 3 re- 
sults of our calculations for Im i;(kn, c^;) at the wave vec- 
tor along the node line of the superconducting order pa- 
rameter in the first BZ. The transition from S(k, uo) oc uo'^ 
to S(k, a;) oc a; for low- lying frequencies is shown for var- 
ious temperatures. Note, the deviation from Landau's 
theory (see solid curve in Fig. 3), Im H ~ cj, results in our 
picture from the strong scattering of the quasi-particles 
on the spin fluctuations and is expected to disappear at 
temperatures T^O, see inset of Fig. 3. In particular, the 
changes in the velocity of quasiparticles are determined 

in EDC as = vp / + versus frequency. At the 

frequencies around 65 meV the Re T^k{^) shows a flatten- 
ing as can be seen via a Kramers-Kronig analysis of Im S. 
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FIG. 4. Calculated frequency dependence of the quasipar- 
ticle self-energy S(ka, a;) at the wavevector k = ka~(l,0.1)7r 
(antinode). The solid curves correspond to the normal state 
at T = 2Tc, whereas the dashed curves refer to the supercon- 
ducting state at T = 0.5Tc. For the wavevector ka the feed- 
back effect of superconductivity on the self-energy is shown. 
Inset: Superconducting gap function A (a;) at wave vector 
k = ka versus frequency. Since the behavior of A and E is 
controlled by Im x(Q5<^) we are able to connect these results 
to the resonance peak observed by INS in cuprates 

Therefore, at this frequency the effect of the scattering 
on spin fluctuations almost disappears. Thus, we find a 
Fermi liquid behavior. Our results also agree with pre- 
vious ones obtained within the spin-fermion model [Q. 
In our microscopic theory we also recover Fermi liquid 
behavior for T~ uj « ujsf- Here, Usf is the characteris- 
tic spin fluctuation energy measured in INS (roughly the 
peak position of Im x(Q, co) |l^]) and is typically around 
25meV for hole-doped superconductors [|2^. Previously 
we have shown that our cogf gives a good description of 
INS data §. On the other hand, for T<Tc the scatter- 
ing is also strongly reduced not only due to Cc; < Usf^ but 
also due to a feedback effect of superconductivity which 
will be discussed in connection with Fig. |^. 

There is a wide discussion whether or not layered 
cuprate superconductors behave like conventional Fermi 
liquids. Earlier experiments (for a review, see Ref. [pT[) 
reveal non- Fermi liquid properties, in particular a lin- 
ear resistivity p{T) for optimal doping, non well-defined 
quasiparticle peaks above the superconducting transi- 
tion temperature Tc seen in ARPES and a strong 
temperature dependence of the uniform spin susceptibil- 
ity observed by nuclear magnetic resonance (NMR) . 
The phenomenological concepts of a Marginal Fermi liq- 
uid (MFL) and a Nested Fermi liquid (NFL) have been 
introduced in order to explain the deviations in the nor- 
mal state from Fermi liquid theory Our results 
shed more light on this question. In agreement with the 
picture of Ruvalds and co-workers we obtain the co- and 
T-dependence of the self-energy mainly due to scatter- 
ing of the quasiparticles on spin fluctuations which is 
strongest for a nested Fermi topology. This also pro- 
vides a microscopic justification for the MFL approach 
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p^ . Thus, for optimal doping {x = 0.15), the micro- 
scopic FLEX approximation includes the phenomenolog- 
ical concepts of both NFL and MFL [|5| . It would be in- 
teresting to extend our studies to the underdoped regime, 
however, the origin of the pseudogap is still unknown [ p6[ . 

In Fig. 4 we demonstrate the feedback of superconduc- 
tivity on S(k, a;). We expect that it is the strongest for 
k?^ (7r,0.l7r) where the gap A{uj) is maximal. One sees 
that mainly the superconducting properties in A(k, cj) 
and in particular in Im A(k, cj) induce changes in the self- 
energy. For the comparison with the experiment we also 
present our results for the superconducting gap. Note, 
this behavior of E(k, cj) and A(k, cj) is related also to 
INS and optical conductivity experiments. In particu- 
lar, the peak position of Im i;(ka, c^) is approximately at 
3Ao — oOsf ^ ^res + ^0 {^res dcuotcs the rcsouaut fre- 
quency observed in INS) according to our previous analy- 
sis ^ . This is in a good agreement with results obtained 
within the frame of the spin-fermion model 

It is remarkable that for electron-doped superconduc- 
tors with a different dispersion Ck |2^, in particular with 
a fiat band lying 300meV below ei? at (7r,0), we get no 
'kink' feature up to frequencies about lOOmeV. This is 
also in agreement with experiment j|] . The reason behind 
this is that Im x(q, ^) has a peak at larger frequencies 
and which is much less pronounced than for hole-doped 
cuprates [ p!9[ . 

In summary, calculating the pronounced momentum 
and frequency dependence of the quasiparticle self-energy 
S in hole-doped high-Tc cuprates we find that this re- 
sults in a 'kink' structure in the dispersion 00 (k) which 
agrees well with recent ARPES experiments. For describ- 
ing the physics in the cuprates it is important that the 
origin of this is the coupling of the quasiparticles to the 
spin fluctuations. The reason for the kink structure is 
a change in the c<;-dependence of the self-energy S from 
non- Fermi liquid to a Fermi liquid behavior. Due to a 
different spectrum Im xi^.^^) of the spin fluctuations in 
electron-doped cuprates we do not find a 'kink' in the 
corresponding spectral density. Furthermore, the feed- 
back effects due to superconductivity on the elementary 
excitations clearly reflect the symmetry of the supercon- 
ducting order parameter. The calculated density of states 
N{uj) = A{uj) = l/NY,^A{kuj) compares weh with SIN 
tunneling data [||. However, due to spatial averaging 
such experiments do not exhibit a kink structure. 
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